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I. INTRODUCTION 

Since the original work of Bender and Boettcher pQ, 
there have been many studies of 'PT-symmetric systems 
and in the last few years the literature on such systems 
has expanded rapidly. In particular, in the seminal pa- 
per mentioned above, the authors showed that a set or 
spectrum of real eigenvalues could be obtained for PT- 
symmetric Hamiltonians of the form H — p 2 + x 2 {ixY, 
the "extended PT oscillator". They found that, when 
e > 0, PT-symmetry is realised as a manifest symmetry 
and consequently all eigenvalues are real; for — 1 < e < 0, 
the PT-symmetry is partially broken, and only a finite 
number of real eigenvalues exist. Below —1, the PT- 
symmetry may be completely broken; however, the situ- 
ation is this case has not yet been examined in detail. 
The feature of a phase transition between a manifest 
and a broken phase is a natural consequence that can be 
expected in any theory displaying two different phases. 
Such a phase transition has been a major focus of in- 
vestigation and details can be found in [5] for this set 
of Hamiltonians. A further discussion can be found in 
Bender's review article, in which several examples are 
presented [3]. Several other papers published since then 
explore the properties of one-dimensional systems, such 
as the PT-symmetric square- well problem [4], or PT- 
symmetric matrix Hamiltonians, see examples in [3]. In 
all of these papers, the aim was to elucidate the point 
that PT-symmetric Hamiltonians, even if they appear 
nonconventional, can display a spectrum of real eigenval- 
ues and that a phase transition from a symmetric to a 
broken phase can occur. 



Recently, it was argued that bosonic field theories de- 
scribed by the (Euclidean) Lagrangian 



c = \m*)) 2 + -^- 2 ^ n+2 w - JbM*), (i) 



where n is a positive integer, can admit several fami- 
lies of real eigenvalues [5]. This raises the question as 
to whether the quantum-mechanical Hamiltonians H for 
the extended PT-symmetric oscillator admit several real 
spectra and not just one for each value of e as shown in [1] , 
and, if so, whether one can classify the number of families 
of solutions that are generated. In this paper, we address 
this question for the Hamiltonians Hi = p 2 + x 2 (ix) e , 
Hn =p 2 + (x 2 ) 5 [6j and Hin = p 2 - (x 2 y , with e, 5 and 
/i > 0, G Z, and thereafter discuss the systematics in the 
field-theoretical case for the Lagrangian in Eq. 0, for 
different values of n. In doing so, we also aim at giving 
a primer on the techniques of complex analysis used, so 
that this paper is self-contained and can be regarded as 
an introduction to this subject. 

This paper is constructed as follows. In Section [TTJ we 
give an overview of the results, review the basic method- 
ology of analysing the asymptotic behaviour of the eigen- 
functions in the complex plane, and discuss how fami- 
lies of real solutions can arise. The analysis is based on 
the WKB approximation [7J. Exact results are presently 
only obtainable numerically and a rigorous proof that 
the eigenvalues are all real is difficult. Analytic contin- 
uation of the eigenvalue problem has been discussed, for 
example, for the anharmonic oscillator in [8]. 

In Section |LTij we apply the methodology to the Hamil- 
tonians Hi, Hn and Hm for positive integral values of 
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e, S and /i, and in Section IV to the generalized field- 
theoretic bosonic Lagrangian, with n taking values from 
to 2 in steps of 1/4. 
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II. BASIC METHODOLOGY 

Symmetry in parity and time is expressed by the 
four relations PxP = —x, PpP = —p, TxT — x and 
TpT = —p, and T is an antilinear operator, TiT = —i. 
Using this prescription, it is easily verified that Hi, Hn 
and Hjjj are all PT-symmetric. It is also obvious that 
the quantum-mechanical Hamiltonians Hjj and Hjjj are 
Hcrmitian in the Dirac sense, i.e. H = Ft. This is also 
true for Hi, for certain values of e. It is the case, for 
example, when e = 0,4,8,...; then Hi e= o = p 2 + x 2 ; 
Hj e =4 — p 2 + x 6 ; Hj e=8 = p 2 + x la , and so on. These 
Hamiltonians all appear bounded from below, and each 
is known to possess a spectrum of real eigenvalues. How- 
ever, by choosing appropriate 'PT-symmetric boundary 
conditions, one can show that these Hamiltonians (with 
e = 0,4,8, ...) in fact have N = (e + 4)/4 different sets 
of eigenvalues or families, which in the WKB approxima- 
tion take on real values. That is, the first Hamiltonian, 
Hi >e —o, has only one real spectrum, Hi te —i may display 
two families of solutions, or have two real spectra, and 
-ff/, e= 8 may have three families with real solutions, or 
three possible spectra of real eigenvalues. 

Continuing to examine Hi, one notes that for e = 
2,6,10,..., these operators take on the form Hi >e —2 — 
p 2 - x 4 , i?/,e=6 = p 2 - x 8 , Hi. t=w = p 2 - x 12 , and 
so on, all of which are Hermitian, but which otherwise 
appear unbounded from below. This set of Hamiltoni- 
ans, however, can also exhibit spectra of real eigenvalues 
when PT-symmetric boundary conditions are imposed. 
Here the number of families that can occur is given by 
N = (e + 2)/4. That is, Hi :e=2 may possess a single 
spectrum of real eigenvalues, iij,e=6 may possess two real 
spectra, while i3j, e =io ma Y have three. 

Finally, the last possible choices of e for Hi, e odd, 
lead to Hamiltonians that are manifestly PT-symmetric 
but not Hermitian: -ffj>=i = p 2 +ix 3 , Hi, e —s = p 2 — ix 5 , 

#1,6=5 = P 2 + ix 7 In these cases N — (e + l)/2 

families appear that can have real eigenvalues. 

The basic methodology underlying these results and 
which is generally applicable to other systems involves 
an analysis of the eigenfunctions in the complex plane. 
Both in PT-symmetric quantum mechanics and in scalar 
field theory, integrals occur whose integration paths can, 
in the most general case, be taken to lie in the complex 
plane. In quantum mechanics, the wave function is 
obtained by integrating the Schrodinger equation and 
imposing boundary conditions, while in field theory, 
the properties of the system can be derived from the 
functional integral formulation of the action in terms 
of the Lagrangian. In each case, one has to analyse 
the convergence properties of the integral; independent 
paths in the complex plane correspond to the possible 
existence of individual families with spectra containing 
real eigenvalues. In the following it is described how, 
with the help of the symmetries of the regions of conver- 
gence of the controlling asymptotic factors, the number 
of real spectra possible both in quantum mechanics and 



in quantum field theory can be deduced. We follow the 
methodology sketched in [3]. 

The first step in dealing with the problem is to iden- 
tify the regions of convergence of the asymptotic wave 
function, known as Stokes' wedges. Figure [l] shows these 
regions associated with the Hamiltonian Hi e=6 = p 2 —x 8 , 
as an example. 




H=p 2 -x 8 



FIG. 1: Regions of convergence for the asymptotic wave func- 
tion associated with Hi, e =6 = p 2 — x 8 (Stokes' wedges). 

As there are always two linearly independent solutions 
of the Schrodinger equation, one finds that the exponen- 
tial behaviour of one solution decays in certain regimes 
(denoted as unshaded Stokes' wedges), while the second 
solution grows in this region. Alternately, the exponen- 
tial part of the second solution decays in the shaded re- 
gions, while the first solution grows. Thus, a wave func- 
tion can always be found whose exponential behaviour 
converges to zero for large \x\ in some wedge of the com- 
plex plane, and the shaded and unshaded regions re- 
flect the difference in sign occurring in the exponent that 
causes the falling or rising asymptotic behaviour. The 
lines separating the wedges are called Stokes' lines and 
these are the lines along which the solution is purely oscil- 
latory. An integration path, however, must lie in regions 
in which the solution converges to zero asymptotically. 
Many paths are possible, but not all necessarily lead to a 
spectrum that can contain real eigenvalues. We conjec- 
ture that the number of paths that can possibly lead to dis- 
tinct real solutions is given by the number of pairs of non- 
contiguous PT -symmetric Stokes' wedges, i.e. Stokes' 
wedges that are symmetric to the imaginary axis, while 
two pairs of Stokes ' wedges which are symmetric to the 
real axis may only be counted once. Hence an integra- 
tion path leads to a spectrum with real eigenvalues only 
if it ends in Stokes' wedges that are PT-symmetrically 
arranged. 

In the case shown in Fig. [TJ following the conjecture 
above, we see by inspection that the system can admit 
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two different real spectra: the shaded pair of Stokes' 
wedges lying under the real axis are T-'T-symmetric to 
one another. This is also the case for the two unshaded 
Stokes' wedges lying directly below the real axis. The 
reflection of this structure into the upper half-plane does 
not bring new solutions. As has already been noted, this 
case is special in that in Hermitian quantum mechanics it 
has no bound states and it only displays a continuous real 
spectrum, as it is not bounded below. As one can see, the 
real axis is not an admissible integration path, reflecting 
this behaviour. In this example, however, we note that 
there are 10 different wedges in the complex plane. In 
the next section we will examine the number of classical 
turning points - in this case there are 8. However, only 
two pairs can be identified that are T-'T-symmetric to one 
another, meaning that two possible spectra of real eigen- 
values can exist. This forms the basis of the conjecture 
made above. 

In scalar quantum field theory the number of real spec- 
tra can be ascertained in a similar fashion. In apposition 
to quantum mechanics, where the sign of the determining 
integral comes from the wave function and its asymptotic 
behaviour always allows for a convergent solution, the 
sign in quantum field theory is fixed. This is due to the 
fact that the determining integral is the Euclidean ac- 
tion, which has a fixed sign. Thus, one finds alternating 
regions of convergence and divergence, limiting the num- 
ber of families of solutions that a bosonic field theory 
can give rise to. Figure [2] shows the pattern of regions of 
convergence and divergence for the Lagrangian of Eq. ([I]) 
when n=2 in one dimension, i.e. for an interaction pro- 
portional to (f> 10 . In this case, the unshaded areas refer 




FIG. 2: Regions of convergence for the n = 2 one dimensional 
scalar field theory (Stokes' wedges), shaded in grey. 



to divergent regions. The integration path may thus only 
lie in the shaded regions. As in quantum mechanics, only 
paths lying in pairs of T-'T-symmetric Stokes' wedges lead 
to real eigenvalues. As a consequence, the number of real 
families of solutions is also given by the number of unique 
pairs of T-'T-symmetric Stokes' wedges. In this case, it is 



seen by inspection that there are three pairs of conver- 
gent wedges that are arranged T-'T-symmetrically. Thus 
three families or spectra can be identified. 

III. FAMILIES OF SOLUTIONS IN 
VT -SYMMETRIC QUANTUM MECHANICS 

In this Section we apply the procedure outlined 
in Section [TT] to the quantum-mechanical Hamilto- 
nians Hi = p 2 + x 2 (ix) e , B. u = p 2 + (x 2 ) 5 and 
Hni = p 2 — {x 2 Y, e, 5, [i > 0, € Z to determine the num- 
ber of possible real spectra or families of real eigenvalues 
that can arise as a consequence of TT-symmetry within 
the WKB approximation [7]. We review the general 
procedure of determining the regions of convergence in 
the complex plane and the energy eigenvalues within 
this approximation. 

As indicated in the last section, our first task is to 
determine the regions in which the wave function con- 
verges in the complex x plane, such as shown for ex- 
ample in Fig. [I] for Hj, e= Q. Determining the regions of 
convergence does not require a detailed solution of the 
associated Schrodinger equation - it is sufficient to inves- 
tigate the asymptotic behaviour of the wavefunction as 
| a; | — > oo, where x is regarded as a complex variable [3]. 
Within the WKB approximation, this limit is determined 
by the exponential controlling factor of the asymptotic 
wavefunction, 

«„~_P( ± /WPW). (2) 

To analyse in which regions of the complex plane this 
function converges to zero when |x| — ¥ oo, one determines 
the lines on which the wavefunction is purely oscillatory, 
i.e. where its real part vanishes. These Stokes' lines are 
thus found by evaluating Eq. |2]) for the appropriate po- 
tential, setting x = rexp(itp), with \x\ — r, and requiring 
that the real part vanish on taking the limit — > oo. 
This condition leads to the Stokes' angle cps; the opening 
angle of each wedge is given by 9 = 2ips- 

To determine thereafter which solution decays or grows 
within each wedge, the anti-Stokes' lines are evaluated 
using Eq. ^ and imposing the condition that the imag- 
inary part of the integral vanishes. This leads to a con- 
dition for ipA- 

Now making further use of the WKB approximation, 
one can determine the real (positive) energy eigenvalues 
associated with each potential, noting that 

Jdxy/E„ -V(x) = (n + l/2)7r, (3) 

X — 

where X- and x+ are the classical turning points evalu- 
ated via the condition E n = V(x). If x is complex, there 
can be several turning points. In the next subsection, we 



4 



will explicitly show that if the complex pairs of turning 
points are chosen to be VT- symmetric, the WKB ap- 
proximation leads to the occurrence of real eigenvalues. 
The number of 'PT-symmetric complex pairs of turning 
points thus determines the number of families containing 
real WKB eigenvalues. 



A. Hamiltonians of the form Hi = p 2 + x 2 (ix) e . 

In the following, we give the results obtained for Hamil- 
tonians of the form Hi = p 2 + x 2 (ix) e . As \x\ — > oo, the 
controlling behaviour of the asymptotic wavefunction is 
given by 



ip ~ exp 



e + 4 

so that the Stokes' lines are given by 



7T 

* - 2 



and the anti-Stokes' lines by 

2k 7T 

y e + 4 2 
The opening angle of the wedges is given by 

6 = 2< / 3 S (fc = 0) 2 



e + 4 



(4) 



(5) 



(6) 



(7) 



The energy spectra associated with the potential 
x 2 (ix) e are obtained by evaluating Eq. To do this, it 
is first necessary to evaluate the turning points and 
x+. Equating E n with the potential (assuming real and 
positive energies), leads to the condition 



E n = -(ix) 



e+2 



(8) 



which has e + 2 possible solutions for x±. In general they 
can be written in the form 



x+ = E% e*>, x.^E'J^e^, 



(9) 



where 7 and £ are angles in which the asymptotic wave 
function converges and which are determined by the po- 
tential. Inserting this ansatz into Eq. (|3| and making 
suitable transformations in the integrals from X— to 
and from to i+, one arrives at the result 



{n+l/2)n = E n t+ ^^ 2t+4 \e^ -e<) / dyy/l - y e + 2 . 

Jo 

(10) 

It follows from this relation that the spectrum can only 
be real if C = n — 7, i-e. if 



_ E l/(e+2) gi(7T-7)_ 



(11) 



Thus the turning points 2;+ and x^ must lie VT- 
symmetrically to one another in order for the eigenvalues 
to be real. In addition, this result is independent of the 
sign of gamma, i.e. it is reflection-symmetric. With this 
choice of C, it follows that 



£(«+4)/(2e+4) 2cos ( 7 ) / (Wl~y £+2 = (n + l/2)7T. 
^0 

Evaluating the integral in the above expression finally 
leads to the closed form 



r(§ 



1 

£+2 



/5F(n+|) 



cos( 7 )r (1 + ^) 



(2e+4)/(e+4) 



. (13) 



with the consequence that the relationship between two 
families of spectra arising from different wedges with con- 
vergence angles 7 and 7' is given by 



cos(7) 
cos (7') 



(26+4)7(6+4) 



(14) 



Examining first the convergence behaviour in the com- 
plex plane, one finds that it differs for the cases in which 
e is odd, even and divisible by four, or even and not divis- 
ible by four. These three cases are illustrated separately. 



The case in which e is odd and the Hamiltonian is 
manifestly non-Hcrmitian is shown in Fig. [3] for e = 
1, 3, 5. When e = 1, and V(x) — ix 3 , there are 

1 /3 

three possible turning points, x± = E n exp(— in/ 6), 

E^/ 3 exp(— j57t/6), E^/ 3 exp(i7r/2). Only one pair 
of these values lies PT-symmetrically to one an- 
other. This pair, cc_ = e}J 3 exp(— i57r/6) and x + = 

1 /3 

E n exp(— in/6), lies within the shaded areas of the 
diagramme, and is thus the only pair that can lead 
to a real spectrum via the WKB calculation. Con- 



tinuing in this fashion, when e = 3 and V{x) = 
—ix 5 , five possible turning points can be identified at 
angles in/10, in/2, -Mir/ 10, -7mt/T0, -IHtt/10. Now 
two pairs are 'PT-symmetric to one another, the pair 
x_ = E\l^ exp(— llwr/TO) and x + = E}J b exp(i7r/10) 
and the pair x_ = E^/ 5 exp(— 7 in/ 10) and x + = 
E n /5 exp(~3i7r/10). The first pair of turning points lies 
in the shaded region of the diagramme in Fig. [3] while the 
second lies in the unshaded pair of wedges below the real 
axis. Each of these pairs of turning points lying in each 
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L-l, H=p*HX } E=3, H=p i -iX s E=S, H=p*+iX ? 



FIG. 3: Regions of convergence (Stokes' wedges) for odd values of e = 1, 3, 5. 



of these pairs of wedges gives rise to a WKB spectrum of 
real eigenvalues. For e = 5, this pattern continues, and 
seven turning points are identified; these can be ordered 
into three pairs that lie 'PT-symmetrically to one another 
and thus lead to three distinct real spectra, according to 
the WKB prescription. 

These results can be summarised into a heuristic pre- 
scription: by examining Fig. [3j we note that for all of 
these cases, one pair of wedges in which the controlling 
factor of the asymptotic solution converges to zero in- 
cludes the real axis and is PT-symmetric, corresponding 
to the generation of one real spectrum. This is the spec- 
trum that was found and discussed in Refs. [T] and [3]. 
In the case when e = 3, another possible real spectrum 
develops, which can be associated with the additional 
pair of 'PT-symmetric wedges (unshaded) that has de- 
veloped in the lower half-plane. For e = 5, yet a fur- 
ther possible real spectrum develops, and this is associ- 



ated with the next pair of noncontiguous PT-symmetric 
wedges (shaded) that has developed in the upper half- 
plane. Thus we note that the number of spectra that 
can form and that is actually determined by the num- 
ber of PT-symmetric pairs of turning points can be in- 
ferred from the convergence pattern of Stokes' wedges in 
the complex plane. The general rule is that this num- 
ber corresponds to the number of noncontiguous VT- 
symmctric pairs of Stokes' wedges. We will simply make 
use of this observation in the next sections. 

To sum up the specific results for odd values of e, we 
find that the general number of real spectra that can be 
generated in this case is given as N = (e + l)/2. Using 
Eq. H we find for e = 1, E n oc (n + 1/2) 6 / 5 ; for e = 3, 
E n oc (n + 1/2) 10 / 7 and the two spectra have the ratio 
R = 1.99 : 1. For e = 5, E n oc (n+ 1/2) 14 / 9 and the ratio 
of the resulting spectra is given as R = 3.52 : 1.41 : 1. 




The case in which e is even and divisible by four, 
i.e. e = 0,4,8,... leads to the Hamiltonians Hj = 
p 2 +x 2 ;p 2 +x 6 andp 2 +£ 10 respectively. The Stokes' and 



anti-Stokes' lines can be easily determined and lead to 
the convergence patterns in the complex plane depicted 
in Fig. [4] In this figure, it is evident that a spectrum 



(> 



of real eigenvalues arises from conventional boundary 
conditions along the real axis, a feature expected from 
Hamiltonians that are Hermitian and bounded from be- 
low. This solution is the only one that can exist when 
e = 0. However, for higher values of e further real spec- 
tra can develop by considering boundary conditions that 
are PT-symmetric. For e = 4 one additional pair of PT- 
symmetric regions develops (from turning points that fall 
within the shaded wedges below the real axis), while for 
e = 8, a pair of shaded wedges and a pair of unshaded 
wedges lie PT-symmetrically to one another, leading to 
two further possible spectra of real eigenvalues. This 
result can be easily verified by evaluating the turning 
points in each case. For e = 4 or V(x) = x 6 , there are 6 
possible turning points, which group into three pairs of 
PT-symmetric values. Since the spectra resulting from 
the turning points in the upper half-plane (in the 'PT- 
symmetric shaded wedges) are the same as those obtained 
from the turning points in the lower half-plane, only two 
distinct spectra are obtained. This reflection symmetry 
in the structure of the upper and lower half-planes is also 
evident for e = 8, for which the potential V(x) = x 10 
yields 5 pairs of PT-symmetric turning points, lying in 
the PT-symmetric Stokes' wedges shown in Fig. [4] Thus, 
in general, the number of distinct spectra containing real 
eigenvalues that may be associated with these Hamilto- 
nians is N = 1 + el A. The lowest spectrum can be calcu- 
lated via Eq. Q to be E n cx n + 1/2, E„ cx (n + 1/2) 3 / 2 
and E n cx (n + 1/2) 5 / 3 , for e = 0,4 and 8 respectively. 
The ratios of the higher order spectra to the lower order 
ones are given as R = 2.83 : 1 and R = 7.09 : 1.43 : 1 for 
e = 4 and 8 respectively. 

Finally, we turn to the case where e is even but not di- 
visible by four, for example, e = 2 or 6. This choice leads 
to the Hamiltonians Hi = p 2 — x A and Hi — p 2 — x s , 
which display the convergence patterns in the complex 
plane that are shown in Fig. [5] and determined from 
Eqs. Q to Although this class of Hamiltonians is 
Hermitian, it is unbounded below. This fact is reflected 
in the figure, as the real axis is a Stokes' line on which 
the controlling factor of the asymptotic wavefunction is 
purely oscillatory: there can be no real (discrete) spec- 
trum arising from boundary conditions given on this axis. 
Nonetheless, there can be spectra containing real eigen- 
values associated with these Hamiltonians: for e = 2, 
the shaded Stokes' wedges represent PT-symmetric re- 
gions in which the controlling factor of the asymptotic 
wavefunction converges to zero, and which can thus give 
rise to a spectrum of real eigenvalues. For e = 6, there 
are two pairs of Stokes' wedges below the real axis, the 
shaded and the unshaded, which lie PT-symmetrically 
to one another and thus lead to two possible spectra of 
real eigenvalues. By symmetry arguments, results from 
the upper half-plane give rise to the same spectra as from 
the lower half-plane. These types of Hamiltonians thus 
have N — (e + 2)/4 possible spectra of real eigenval- 
ues. For the cases displayed in the figure, the discrete 
energy spectrum for the Hamiltonian e = 2 is given as 



E n cx (n + 1/2) 4 / 3 , while for e = 6, E n cx (n + 1/2) 8 / 5 
and the relation of the higher order solutions to the lower 
is R= 4.10 : 1. 





i-6, H=p*-X* 



FIG. 5: Regions of convergence (Stokes' wedges) for even val- 
ues of e that are not divisible by four, for example e = 2, 6. 



B. Hamiltonians of the form Hh = p 2 + (x 2 ) 6 . 

Hamiltonians discussed in the previous section fall for- 
mally into one single class, Hi = p 2 + x 2 (ix) e . Viewing 
them individually for even values of e, one notes that the 
sign of the potential alternates, a feature that gives rise 
to a rich behaviour. For real values of the potential, it is 
interesting to note how a change of sign alters the distri- 
bution of possible real (discrete) spectra. For this reason, 
we evaluate the spectra for Hamiltonians that are deter- 
mined by Ha = p 2 + (x 2 ) s , with a view to examining 
those Hamiltonians that are not described by Hi for any 
value of e. Now the controlling factor in the asymptotic 
form of the wavefunction from Eq. ^ behaves as 



ij){x) ~ exp I ± 



„<5+l 



1 



(15) 



in the large x limit. The Stokes' lines then follow on 
regarding x as a complex variable and setting the real 
part of the wavefunction to zero. One finds 



6 + V 



(16) 



where k € Z. The anti-Stokes' lines occur at angles 



1 



(17) 



We illustrate the results for the cases 5 = 2 and 4 in 
Fig. [6j which differ markedly from that of the previ- 
ous Hamiltonians of Subsection A: in both cases PT- 
symmetric Stokes' wedges lie along the real axis. Each 
wedge, however, arises from a different convergent so- 
lution and is therefore differently shaded. The combi- 
nation, as is well-known in quantum mechanics, leads 
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to a spectrum of real eigenvalues. For 6 = 2 and 
Hi i = p 2 + x 4 , this is the only real spectrum that exists. 
For 8 = 4, a second pair of PT-symmetric wedges has 
developed in which the wave function converges asymp- 
totically to zero, leading to a second possible spectrum 
of real eigenvalues for Hn = p 2 + x 8 . In general, the 
number of families found for even 8 is N = 8/2. The 
case 8 odd, for example 5 = 1,3, maps onto the cases 
already discussed in Section A, with e = 0, 4. The num- 
ber of spectra of real eigenvalues obtained in this case is 
N = (8+l)/A. 




b-3, H=p 2 +X 1 6=4, H=p'+X> 



E n w 2.81(71 + 1/2) 3 / 2 . For odd values of fx, the number 
of real spectra that exist is given by N = (/x— l)/2. Even 
values of \l correspond to the cases e = 2 and 6 respec- 
tively of Hi detailed in Section A, leading to the result 
that the number of real discrete spectra in this case is 
N = n/2. 



Im k 



.21 




11=1, H=p-'-x l n=3, H=p l -x f 



FIG. 7: Regions of convergence (Stokes' wedges) for (i = l,3, 
representing the Hamiltonians H = p 2 — x 2 and H = p 2 — x 6 . 



FIG. 6: Regions of convergence (Stokes' wedges) for 8 
representing the Hamiltonians H 



p 2 + x 4 and H 



2,4, 



p 2 +x° 



C. Hamiltonians of the form Hiii — p 2 — (x 2 ) 11 . 

One distinct case in the type of Hamiltonians that can 
be constructed and which has not yet been dealt with 
fully is given by Hm = p 2 — (x 2 )^ . Here, the controlling 
behaviour of the asymptotic wavefunction is given by 



ip(x) ~ exp ( ±i 



1 



II- 



1 



(18) 



for large x. Regarding x as complex leads to the Stokes' 
lines occurring at angles 



fs = 



k'TT 
fl+1 



and the anti-Stokes' lines at 



k- 1 

2 



/i+ 1 



(19) 



(20) 



The cases fi — 1 and /i = 3 are illustrated in Fig. [7] In 
both of these cases, as with all the quantum-mechanical 
Hamiltonians, regions in which the controlling factor of 
the asymptotic wave function converges to zero exist in 
the entire complex x-plane. However for fi = 1, there are 
no noncontiguous 'FT-symmetric wedges. In this case, no 
spectra with real eigenvalues exist. For the case /i = 3, 
one pair of noncontiguous "PT-symmetric wedges devel- 
ops lying just below the real axis (left shaded, right un- 
shaded). Here one spectrum with real eigenvalues can 
arise. In this case, the WKB approximation leads to 



IV. FAMILIES OF SOLUTIONS IN 
VT -SYMMETRIC SCALAR FIELD THEORY 

In this Section we review the analysis of 'PT-symmetric 
scalar quantum field theory, following [5]. That is, we 
show how to obtain the regions of convergence of the 
partition function in the complex plane, explain how to 
deduce the number of real spectra from this image and fi- 
nally obtain the particle masses that are associated with 
each spectrum. In doing so, it is our aim to clarify the dif- 
ference between the quantum-mechanical cases described 
in the previous section and the field-theoretical situation. 

We investigate systems with the set of Lagrangians 
given by 



C 



1 



(d^x)) 2 



fJ 



An + 2 



<p in+2 {x) - J(x)(j)(x) 



for different values of n. The partition function associ- 
ated with this Lagrangian is given by 

Z[J(x)} = (0|0) =J c V ( />(x)exp{- J d D x'[\{d<t>(x')) 2 

+ 4^^"+ V) - W)]}, 

where D is the dimensionality. The simplest case which is 
in direct analogy to the quantum-mechanical one is that 
in one dimension. The controlling behaviour of Z[J(x)] 
is then determined by the integral 



^exp{-(^ 4,l + 2 } 



(21) 



In a similar fashion to quantum mechanics, one can 
analyse for which values of ip the integral converges to 
zero when \ip\ — > oo. 




FIG. 8: Regions of convergence (Stokes' wedges) for n = 0, 1, 2 



Furthermore, it is possible to determine the energy 
spectra associated with C and thus the particle masses. 
To do so, we follow the treatment given in [5]. The field 
equation is deduced from the Lagrangian to be 

~d 2 (/){x)+g^ n+1 {x) = J{x). 

Taking the expectation value of this equation and using 
the defining relations for the vacuum expectation values 
and the connected Green functions, 

- 5 n Z[J{x)] 6"\nZ{J(x)} 
W W) - 5Jn{x) and G n{x,x, ...) - 5jn{x) , 

one can derive the infinite series of coupled Schwinger- 
Dyson equations. Truncating Green functions of third 
and higher order, one finally obtains two equations for 
Gi and G2, the latter of which determines the mass M: 



He 



"2n+l 



= 0, 



(22) 



where Hem+\ are Hermite polynomials of order In + 1, 
and 



(-V 2 +M 2 )G 2 (x-y) = 6 D (x-y), (23) 



with 



M 2 = (G 2 (0)) 2n He' An+l (id/V^O)) . (24) 



Solving Eqs. (22) - (24 1 for M and G\, one generally 



obtains more than one real and positive solution for 
the mass M, while G\ can be 0, or negative and purely 
imaginary. This latter case represents the "PT-symmetric 
solution of the system. 

For completeness, we consider first the case when n 
takes on integral values, i.e. n = 1,2,3.... The region 
in which Eq. (21) converges is again determined by the 



Writing if — |</?|e 41? , it follows that these angles occur at 




and 



,'k + l/2 s 
1 1„ • 2 1 



4n 



(25) 



(26) 



Stokes' or anti-Stokes' lines associated with this integral. 



where k is an integer. The controlling parts of the field- 
theoretical integrals - unlike their quantum-mechanical 
counterparts - do not converge to zero in all regions of the 
complex (p plane. The convergence patterns are shown in 
Fig. [8] for the cases n — 0, 1, 2. Each case has a conver- 
gent region along the real axis, yielding one real energy 
spectrum, or, for the ground state, the particle mass. In 
addition, however, the case n = 1 develops further re- 
gions of convergence, which lie ■pT-symmetrically to one 
another and thus generate a further spectrum leading to 
a further mass that has the ratio R — 1.62 : 1 to the 
first mass. When n = 2, a second extra real energy spec- 
trum develops, so that there are three in total, leading 
to masses lying in the ratio R — 2.26 : 1.19 : 1 to one 
another. The region of convergence for n — strongly re- 
sembles the quantum-mechanical case depicted in Fig. [4] 
when e = 0. The fact that there is no convergent solu- 
tion in the field-theoretic case in the vertical wedges does 
not play any role in determining the number of real spec- 
tra. These two cases are equivalent. For the case that 
n = 1, however, the distribution of Stokes' wedges shown 
in Fig. [8] resembles the case e = 8 in Fig. |3J but has one 
less available pair of PT-symmetric wedges to construct 
a real spectrum. The quantum-mechanical case thus ad- 
mits one more real spectrum than the field-theoretical 
case with a similar wedge structure. The number of real 
spectra in the field-theoretical case is given by N = n + 1. 

For half-integer values of n, such as n — 1/2 or 
n = 3/2, the picture shown in Fig. [9] emerges. For 



9 



n = 1/2, only one real energy spectrum can exist, while 
for n — 3/2 an additional "PT-symmetrical spectrum can 
be found, with mass ratio R — 1.31 : 1 between the two 
solutions. In general, for half-integer values of n there 
are N = n+ 1/2 possible real spectra. 




n=l/2; tp- interaction n=3/2; <p s interaction 



FIG. 9: Regions of convergence (Stokes' wedges) for n = 1/2 
and 3/2. 

Finally, we note that there are no PT-symmetric so- 
lutions for values of n = 1/4, 3/4, 5/4, . . . corresponding 
to interactions of the form (p 3 , ip 5 and ip 7 respectively. 
This can be seen immediately on examining the pattern 
of Stokes' wedges for example, for n = 1/4, shown in 
Flo- [Tn] 





Im 














Re 













FIG. 10: Regions of convergence (Stokes' wedges) for n = 1/4. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have described the classification of 
families of (positive) real spectra associated with the 
quantum-mechanical Hamiltonians Hj — p 2 + x 2 (ixY, 
H u = p 2 + (x 2 ) s and H UI = p 2 - (x 2 )'*, with e,S and 
(i > 0; e, 5, }i £ Z within the WKB approximation and for 



a scalar field theory £ = \{d(j){x)) 2 +g/{An+2)(j) in+2 {x)~ 
J(x)(f>(x) within the functional integral formalism. In 
the quantum-mechanical case, we have shown that the 
number of real spectra depends on the number of pairs 
of PT-symmetric turning points that can be found in 
the complex plane. This in turn depends on the type 
of Hamiltonian and the form of the interaction poten- 
tial. The classification of the number of PT-symmetric 
turning points can be reinterpreted in terms of the posi- 
tion and number of Stokes' wedges in the complex plane, 
leading to the rule that each pair of noncontiguous PT- 
symmctric Stokes' wedges may result in a positive energy 
spectrum containing real eigenvalues; two pairs of Stokes' 
wedges which are symmetric to the real axis may how- 
ever only be counted once. Thus, by simply examining 
the Stokes' wedges in the complex plane, one can deduce 
the number of spectra of real eigenvalues that can oc- 
cur. One recovers the usual result that Dirac-Hermitian 
Hamiltonians have a spectrum of real eigenvalues due to 
"normal" boundary conditions placed along the real axis. 
However, these Hamiltonians can give rise to additional 
real spectra if additional PT-symmetrical regions of con- 
vergence exist in the complex plane. Non-Hermitian, but 
PT-symmetric Hamiltonians can also possess a spectrum 
of real eigenvalues (in fact several), while the class of 
Dirac-Hermitian Hamiltonians that appear unbounded 
from below can also produce several such spectra or fam- 
ilies. This is because the Hamiltonians must always be 
solved together with a set of given boundary conditions. 
For all cases, the spectra can be evaluated simply within 
the WKB approximation and the ratios of the different 
spectra can be found. 

A similar analysis of the scalar field theory using func- 
tional integrals shows that again the structure of the 
Stokes' wedges in the complex plane suffices to give infor- 
mation on the number of possible real spectra. However, 
fewer real spectra are available in cases which appear for- 
mally similar to their quantum-mechanical counterparts, 
as the regions of convergence of the asymptotic control- 
ling factor in the complex plane are restricted by the sign 
of the power of the action, as it occurs in the generating 
functional. The ratio of the masses (ground state ener- 
gies) of the individual spectra can also be simply eval- 
uated. Here the approximation lies in the truncation of 
the Dyson-Schwinger equations. 

Finally, it is appropriate to stress the limitations of the 
analysis given here. In the quantum-mechanical cases, 
it is unclear how well the WKB approximation repro- 
duces the actual energy eigenvalues, which can only be 
calculated definitively numerically. Usually the result is 
excellent for the lowest states, see [3]. It is, however, 
not possible to use it to describe the spontaneous break- 
down of PT-symmetry, and the associated phase transi- 
tion. Nor can it be used to address deeper questions such 
as whether a particular spectrum has a finite or infinite 
number of real eigenvalues. More powerful techniques are 
needed. Some useful references that go beyond the WKB 
approach can be found for example in and |10j . 
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